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Abstract 

A higher dimensional modified gravity theory with an action that 
includes dimensionally continued Euler-Poincare forms up to second 
order in curvatures is considered. The variational field equations are 
derived. Matter in the universe at large scales is modeled by a fluid 
satisfying an equation of state with dimensional dichotomy. We study 
solutions that describe higher dimensional steady state cosmologies 
with constant volume for which the three dimensional external space 
is expanding at an accelerated rate while the (compact) internal space 
is contracting. We showed that the second order Euler-Poincare term 
in the constructions of higher dimensional steady state cosmologies 
could be crucial. 
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1 Introduction 


The idea of considering an early inflationary epoch during the evolution of 
the Universe that is characterized by an exponential expansion (that lasted 
till 10 35 seconds after the Big Bang, corresponding to energy scales 
~ 10 16 GeV) is successful in resolving the problems of standard Big Bang 
cosmology such as the horizon and flatness problems. Furthermore the same 
idea provides strong clues for an explanation of the origin of large scale 
temperature fluctuations observed in the cosmic microwave background ra¬ 
diation p la in g. The actual models of inflation have several variants that 
are mostly based on general relativity in which the inflation would be driven 
by scalar field(s) with ad hoc potent.ial(s); see for instance Ref.[5] for a recent 
review. We do not yet have a concrete and unique realization of inflation from 
a fundamental theory. Leaving aside the accelerated expansion in the early 
stages of the universe, it is today confirmed beyond any doubt by means 
of several independent observations that the universe has again started to 
expand at an accelerated rate approximately 6 x 10 9 years ago. A satisfac¬ 
tory explanation of this current acceleration that takes place at energy scales 
lCU 4 eV where we supposedly know physics very well is lacking. The most 
successful cosmological model accommodating the observed pattern of ex¬ 
pansion of the universe so far is the six parameter base ACDM model that is 
simple and in reasonably good agreement with available high precision data 
m®. Yet two well-known problems related with the A assumption, namely 
the coincidence and fine-tuning problems, may be signaling that it should 
rather be regarded as a good approximation or a limiting case of a more 
general model in which the vacuum energy need not necessarily be realized 
in terms of a cosmological constant A and may deviate from A considerably 
in the far past and/or future [HI El HO, Hj. This calls for a dynamical de¬ 
scription of dark energy and/or a possible modification of general relativity. 
Some recently reported tensions between the A assumption and high preci¬ 
sion data could also be resolved in case of evolving dark energy (DE) and/or 
a modified theory of gravity EU El Eg. 

It is of course more desirable if inflation, or evolution scenarios that can be 
alternative to inflation, and the late time acceleration of the universe can 
all be derived from first principles or from a unified theory of fundamental 
interactions of nature such as the string theories or M-theory. A common 
property of such unified theories is that their consistent constructions re- 
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quire the existence of extra spatial dimensions in addition to our familiar 
four space-time dimensions. It is generally argued that we live in a universe 
in which all but four of the space-time dimensions are compactified on an 
unobservable, non-singular internal space without boundary; thus leaving 
behind an observable (3 + l)-dimensional external space-time. But it is also 
known that the accelerated expansion of our external space is difficult, to 
maintain in a cosmological context with such an approach unless the dimen¬ 
sional reduction of the internal space is dynamical mmm- It should not 
be surprising that considerations of extra dimensions also led to generaliza¬ 
tions of Einstein’s gravity in a natural way [191 12Q1[21] [221 [23] [2U [25J 1261127] . 
In fact long before the string models, specific higher dimensional generaliza¬ 
tions of the Einstein tensor were pointed out such that they are covariantly 
constant and contain at most second derivatives of the metric components 
[28, [29] . These tensor concomitants follow by a variational principle from 
the dimensionally continued Euler-Poincare densities. The standard Kaluza- 
Klein reduction of gravitational actions that include dimensionally continued 
Euler-Poiucare forms exhibits a complicated pattern of non-minimal inter¬ 
actions in four dimensions [3U1EU E2]- Their various cosmological solutions 
have also been much studied [33] [34, 35j [36, 37J [38]. However, the relatively 
recent notion of dark energy invoked a renewed upsurge of interest in their 
cosmological solutions [39] 00] 0H 021 03] 0H 05]. 

In this paper, for the sake of simplicity, we are going to neglect the dilaton 
and axion fields that appear in effective string field theories. On the other 
hand, a cosmological const ant/vacuum energy is not expected to occur in the 
effective low-energy theories of superstrings, but its existence is not excluded 
and it typically predicts negative vacuum energy, say a negative cosmological 
constant. Accordingly we consider also a negative cosmological constant 
in our model since its presence will lead to more general models without 
complicating the field equations much. We are particularly interested in the 
dynamics of our external space in a higher dimensional steady state universe 
that was characterized by two main assumptions in [46] 147] : (i) the higher 
dimensional universe has a constant volume as a whole but both the internal 
and external spaces are dynamical, (ii) The energy density is constant in the 
higher dimensional universe. In this paper, we consider a higher dimensional 
fluid to model the matter distribution in the universe that exerts pressures 
that are distinct along the external and internal dimensions in accordance 
with the first assumption. Introduction of a higher dimensional anisotropic 
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fluid provides us with new degrees of freedom that would in return allow us to 
deal with the dimensional dichotomy of our higher dimensional steady state 
universe without over-determining the system of field equations. We will be 
relaxing the second assumption accordingly and instead study a more general 
case in which the equation of state parameters of the higher dimensional fluid 
along the external and internal dimensions are dynamical but the difference 
between them is kept constant. 

In section ([2]), we give the gravitational field equations derived by a vari¬ 
ational principle from a truncated action that includes only up to second 
order Euler-Poincare forms in the gravitational sector. The vanishing of the 
space-time torsion is handled by considering Palatini-type variations using 
the method of Lagrange multipliers. In section (J2D the cosmological space- 
time geometry is specified through a spatially homogeneous and flat but not 
necessarily isotropic (1 + 3 + n)-dimensional synchronous space-time metric 
in locally Cartesian coordinates. We also write explicitly the stress-energy- 
momentum tensor of a higher dimensional fluid. The reduced field equations 
are found. In section (0]), the constant volume condition is imposed and the 
field equations are further reduced. We also discuss an equation of state that 
respects the dimensional dichotomy of the spatial geometry. Then exact solu¬ 
tions are found and classified according to the ranges of our free parameters. 
Section (J5]) is devoted to concluding remarks. 

2 Variational Field Equations 

Einstein’s gravitational field equations in 4-dimensions is a non-linear system 
of second order partial differential equations that can be derived by a varia¬ 
tional principle from the Einstein-Hilbert action that involves a Lagrangian 
density linear in the curvature components. It is not equally well-known 
that in dimensions higher than four, a gravitational Lagrangian density can 
be supplemented by unique additional terms, besides the Einstein-Hilbert 
term and a possible cosmological constant, without destroying the desirable 
property that the variational field equations do not carry derivatives of the 
metric components of order greater than two. Such additional terms in the 
action are essentially topological in nature and are given by the integrals of 
dimensionally continued Euler-Poincare forms over the space-time manifold 
M. The zeroth order Euler-Poincare form is taken to be the volume form it¬ 
self and yields the cosmological constant. The first order Euler-Poincare form 
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is linear in the curvature components and is equal to the Einstein-Hilbert La- 
grangian density. In D = 2 dimensions this is a closed form whose integral 
over a 2-manifold M with boundary leads to the famous Gauss-Bonnet the¬ 
orem. In dimensions D > 2 its variations lead to Einstein held equations. 
In a similar way, the second order Euler-Poincare form that is quadratic in 
the curvature components is a closed form in D = 4 dimensions. Its integral 
over a 4-manifold M gives the topological Euler number of M. The second 
order Euler-Poincare density when continued to lower dimensions (D < 4) 
vanishes identically, while when continued to higher dimensions (D > 4) its 
variations lead to second order dynamical held equations. The higher dimen¬ 
sional modified gravitational action obtained by integrating the second order 
Euler-Poincare form over M is sometimes called the Gauss-Bonnet gravity in 
current literature. It was further pointed out by Lovelock long time ago that 
the n-th order Euler-Poincare form is a closed form in D = 2n dimensions 
and if used as part of a gravitational action in dimensions D > 2n contributes 
only in second order to the highly non-linear variational held equations. 

Thus the most general gravitational Lagrangian density in D -dimensions 
that yields second order variational held equations can be written as a linear 
combination of all possible dimensionally continued Euler-Poincare densities: 


D — l 
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( 1 ) 


where we introduced (dimension-full) coupling constants A n and the Euler- 
Poincare forms are given by 



( 2 ) 


R a b = duj a b + 1 o a c A a f b are the Riemann curvature 2-forms of (torsion-free) 
Levi-Civita connection 1-forms {u“J. They are uniquely determined by a 
Lorentzian signatured metric tensor g = r] a be a <g)e b where r] a b = diag(— +++), 
given in terms of an orthonormal co-frame {e a }, through the Cartan structure 
equations de a + u> a b A e b = 0. We further introduced the notation e ab '" c = 
e a A e b A ... A e c and the Hodge *-map that hxes the space-time orientation 
with the volume H-fornr *1 = e m2 --- D ~ l _ 

Here we vary the total action (in D > 5 dimensions) 



( 3 ) 
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where we take a truncated gravitational Lagrangian density given by 

L g = A 0 * 1 + l -R ab A *e ab + ^ R ab A R cd A *e abcd , (4) 

where Ao stands as a higher dimensional negative cosmological constant. The 
constraint Lagrangian 


L c = [de a + u a b A e b ) A [i a (5) 

imposes the zero-torsion constraint on the connection through the variations 
of the Lagrange multiplier (D — 2)-forms /i a . k > 0 denotes the Newton’s 
constant that is scaled appropriately in 74-dimensions. Then the co-frame 
variations of the action (|3]) give 

A 0 * e a + ^R bc A *e abc + -jR bc A R df A *e abcdf + Dn a + KT a (m) = 0, (6) 

while the connection variations of (J3]) give 

*e abc A T c + y * e abcdf A R cd A T f + KY> ab (m) = e a A fi b - e b A /i a . (7) 

In these expressions we introduced through the variations of the matter La¬ 
grangian, the stress-energy-momentum tensor T ab {m) and the angular mo¬ 
mentum tensor S abtC (m ) to be read from 

(“5“?) = r “( m ) = T ab{m) * e b ( 8 ) 

and 

= = Sab ' c ( m ) * ( 9 ) 

respectively. The variations of the action with respect to the Lagrange mul¬ 
tipliers impose the zero-torsion constraint T a = 0 and the held equations 
OS) and (J7j) above are to be solved subject to this constraint. Then the 
connection-variation equations simplify and an algebraic solution for the La¬ 
grange multiplier forms can be found as follows: 

ha = K (^X^a ~ 2^ A tXciX b^) ' ( 10 ) 
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We use interior products such that Av 0 (e 6 ) = b b a with respect to orthonormal 
frame vectors {X a }. Then the Lagrange multipliers are substituted into the 
remaining field equations and we arrive at the final form of our modified 
Einstein field equations 

Ao * e a + -R bc A *e a bc + ~^R bc A R df A *e abcd f 

= - KT a (m ) + nD{b Xi ^a) + ^e a A D{b Xb b X Jl bc ). (H) 

As a last remark we note that the integrability of the Einstein field equations 
(fTTj) yields the conservation law of matter: 

Dr a (m ) = R ab A b Xc E cb (m). (12) 


3 Cosmological Ansatz 


The geometry of a cosmological space-time will be given by the metric 

3 n 

g = -dt 2 + R{t) 2 ^(dx 1 ) 2 + S(t) 2 ^2(dy a ) 2 (13) 

i =1 a=l 

where {x 1 } are the Cartesian coordinates of a flat (external) 3-space and {y a } 
are the coordinates of a (locally) flat, compact (internal) n-space. R(t) and 
S(t) are the corresponding scale factors of cosmic time t. Let us choose the 
following orthonormal basis 1-forms 

e° = dt, e* = R{t)dx l , e a = S(t)dy a , (14) 


so that the non-zero connection 1-forms are determined to be 

o j q S 

<jy = — e <jy = — e 

W 1 R ^ ! ^ a S ^ ’ 


(15) 


where • = 4- denotes time derivative. Therefore the curvature 2-forms are 

at 

given by 


R R 


S S 2 


R\ = (4 + id eU Ae\ R\= [4 + ^)e°Ae a , 


R R 2 


S S 2 


R 


S 2 


RS 


R l j = ^e*Ae», R a p = ^e a AeP, R\ = -- e 'Ae a . ( 16 ) 
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We will model the matter distribution in a D- dimensional space-time by a 
homogeneous but anisotropic fluid with no rotation. Thus we take E a6 (m) = 
0 and the stress-energy-momentum tensor of the fluid will be specified by 

3 n 

T = T ab e a <g> e b = p(t)dt 2 + p ex t(t) dx l ) 2 + p int (t) dy a ) 2 (17) 

2—1 a =1 

where we assume that the mass density p > 0, but the pressures p ex t(t ) and 
Pi n t(t ) are arbitrary except for the conservation law they should satisfy: 

R S 

P + 3-^{p + Pext)+n-(p + p int ) = 0. (18) 


At this point, we skip calculational details and give the independent equations 
of motion we obtained as follows: 


Aq + 3 


R 2 0 RS n(n — 1) S' 2 1 
R? + RS + 2 S 2 J 


+ A 2 |9n(n-1) — — 

d3 C d 03 

+ 6 „_- + 3n(n-l)(n-2)-- 

77,(77. — l)(n — 2) (77 — 3) S' 4 1 
+ -4-Si) 


(19) 


Ao 


n R R 2 S n RS 77(77 — 1) S 2 

2 R + m +n s + 2n Rs + ^—& 


R 


s 2 


RS 


A 2 y 2 n(n - 1)— + 4n— - 


R 


S 2 


RS 


S 


R 


RS 


S 2 


- 2n— + 4n(n - 1) —— + n(n - l)(n - 2) — 


S 


R 2 


RS 


S 2 


732 c2 73 03 

+3n(n - 1)^ + 2n(n - 1 )(n - 2)|^ 
77(77 — l)(n — 2) (77 — 3) S' 4 1 

^ g4 r = ~ K Pext j 


( 20 ) 
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( 21 ) 


The coupled equations (fT9l) . (T20l) and (|2B are to be satisfied by five functions 
R , S, p, p ext and pi nt and therefore the system is not fully determined as it 
is. Two extra constraints must be provided by some extra assumptions in 
order to get fully determined solutions. 

4 Exact Solutions with Constant Volume 

We are particularly interested in higher dimensional steady state cosmologies 
characterized by the following two properties: (i) the higher dimensional uni¬ 
verse has a constant volume as a whole but the internal and external spaces 
are dynamical, (ii) The mass density is constant in the higher dimensional 
universe. These are strong assumptions that had been motivated and dis¬ 
cussed in the context of similar cosmological models in our previous works 
pi sn. We will further add comments towards their justification in the 
concluding section E3 Here it should be noted, however, that we allow the 
energy density to be dynamical in a certain way, so that a constant mass 
density can be studied as a special case. Therefore, firstly we assume that 
the (3 + n)-dimensional volume scale factor of the universe is constant: 


R 3 S ” = v 0 . 


( 22 ) 


It will be convenient to introduce the Hubble functions of external and inter¬ 


nal spaces given by H = ^ and H' — |, respectively. Then it follows from 
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(122]) that 3 H = —nH', assuring the dynamical contraction (hence reduction) 
of the internal space (H 1 < 0) for an expanding external space (H > 0). The 
held equations in terms of H read: 

3t7 —I— Q 

A 0 - H 2 + 3\ 2 b n H 4 = K p, (23) 

2 n 


• 3 71 -1 - Q 

Ao — H -\ -—— H 2 + \ 2 H 2 ( a n H — b n H 2 ) = —np ext: 

3 rV 3n + 9 „ 2 


An + -H + 


n 


2 n 


H 2 + A 2 H 2 (a n - 4 b n )H - b n H 2 = - K p int 


where we se10 


9 (n 2 + 3n — 6) ; 3(n + 3)(n 2 + 15n — 18) 

^n. “ • 0 n 


n 


4n 3 


Furthermore our conservation law of matter reduces to 


(24) 

(25) 


(26) 


P “ 1 “ 3 H (Pext Pint ) 0 * 


(27) 


Secondly we assume a simple equation of state 


Pext Pint _ 

P 


(28) 


where Aw is a constant that we propose to call dimensional dichotomy pa¬ 
rameter. It follows that p oc Rr :iAw in general. A higher dimensional isotropic 
distribution of matter with constant mass density follows in the case Aw = 0. 
If there is a dimensional dichotomy i.e. Aw ^ 0, on the other hand, then, 
as the external space expands, the mass density decreases if p ext > Pint and 
increases if p ext < p in t- 

In order to solve this system we proceed in the following way. We differentiate 
both sides of (|23|) and on the RHS use the conservation relation (1271) together 
with our equation of state (1281) to replace p. Then we use the first equation 
(1271) again to eliminate p itself, arriving at a first order master equation 
satisfied by H. We will solve this master equation and classify its solutions 
according to the signs and magnitudes of the cosmological coupling constants 

1 Note that a n and b n are both negative for n = 1 and both positive for all n > 2. For 

the sake of simplicity we will keep the second case in what follows. 
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Ao and A 2 . Once an expression for H is fixed, the external and internal 
pressure functions can be found by substituting this H into the equations 
(|24|) and (|25|) . respectively. Finally we integrate H and determine the scale 
factors R and S, thus completely specifying the cosmology. 

From now on we will be looking at the master equation for H that reads 


3n + 9 A 
2 n ) 


+ 


n[n 


+ 3) 


'3n + 9 


2(n 2 + 15n-18)A 2 J \ 2 n 


H 



.4 n(n + 3) 1 

~ (n 2 + 15n - 18) W 


1 2 4- 15 n - 181 A 


/3n + 9A 

V 2n ) 


H 2 

(29) 


We re-write it as a simple differential equation 


x = OJ 



A) 


Q 

(:x 2 - A). 


(30) 


where we dehned 


'3 n + 9 

x=\i - H, uj = 

2 n 


2 n 3Aw 
3n + 9 4 


(31) 


with two free parameters 


^ _ n[n + 3) 


(n 2 + 15n- 18) 2A 2 ’ 


Q = \ 2 - 2A 0 A. 


(32) 


It is worth to note that the held equation (1251) with the above re-parametrization 
gives for the energy density 


^ [ 4 2 - A ) 2 - Q] ■ 


(33) 


We first consider briefly the case of an isotropic distribution of matter. Sup¬ 
pose A 7 ^ 0. We take uj = 0 and the master equation becomes x = 0. This 
leads to an exponentially expanding external scale factor R(t) = R(0)e Hot 
and an exponentially contracting internal scale factor S(t) = S(0)e~™ Hot . 
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For isotropic solutions with A ^ 0, the energy density may be positive def¬ 
inite or negative definite or zero depending on the roots of the algebraic 
equation 


3n + 9 2 1 /3n + 9\ 

~^r Ho + 2x{~^r) 


Ht 


0 . 


Next we consider cases with uj ^ 0 and A ^ 0: 

(1) Suppose Q ^ 0. The generic solution is implicitly given by 


Va + VQ - x 

1 

2y/X+VQ 

1 

I 

> 

1 

Wx-VQ 

y/ A + y/Q + X 


y/ A — y/Q + x 



(34) 


where C is an integration constant. We cannot further develop such solutions 
and obtain explicit expression for R(t) and S(t) in general. 

( 2 ) Suppose <3 = 0 , that is A = 2A 0 . Now the generic solution reads 


Vx- 


x 


Vx + 


x 


i 

Vx 


= Ce 


2 ut 


(35) 


and can be inverted to write down explicit expressions. For A > 0 we have 


f — VXcoth(\/Xujt), x 2 > A 
\ — y/\ tanh(\/Auk), x 2 < A 


(36) 


These expressions can be integrated for the scale factor of the external space 
so that 

f | sinh(vXnt)| 4 /( 3Au, \ x 2 > A 
\ | cosh(\/Aa;t)| 4 ^ 3Au '\ x 2 < A 

The corresponding mass density function turns out to be positive definite: 




f | skill 4 (a/Ac ot), x 2 > A 
\ | cosh~ 4 (\/Acut), x 2 < A 


(38) 


For A < 0 we have 

x(t) = — \/|ATcot(-\/|A[a;t). (39) 

Again we integrate the above expression and find that 

R(t) = | sin(vlAjwt) | 4/(3Au,) . (40) 
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The corresponding mass density function for this solution is negative definite: 


np(t) = —-—-sin 4 (\/Xu}t). 


(41) 


We also wish to point out some special classes of solutions which cannot be 
obtained from above as limiting cases. 

(A) First consider the case of Einstein-Hilbert action with a cosmological 
constant where the Euler-Poincare term is absent. We set A 2 = 0 accordingly 
but still keep oj 7 ^ 0. Then the master equation reads 


x = 2 uj ( x 2 — A 0 ) 

For Aq > 0 the solutions are given by 


x 


(t) = 


— a/Ao coth(2-\/Aou;£), x 2 > Ao 

—v^Aq tanh(2-\/Aoa;t), x 2 < A 0 


The external scale factor can be calculated as 

cosh(2^X~ 0 ujt)\ 2 ^ Aw \ x 2 > A 0 
sinh( 2 ^a;f)| 2 /( 3A -), x 2 < A 0 


m = 


which give a positive definite the energy density for 0 < A 0 < x 5 

Kp{t) = - / _ _ 7 ) 

cosh 2 (2y/\^u>t) 

and a negative definite energy density for x 2 < A 0 

Ao 


np{t) = 


sinh 2 (2i/A(Wt) 


For An < 0 we have 


x(t) = y/ / |Ao| tan(2\/|Ao[a;t) 


and 


R(t) = |cos( 2 v ^cuf)| 2/(3A!i ’ ) . 

In this case, the corresponding energy density is always negative definite: 

l^o | 


(42) 

(43) 

(44) 

(45) 

(46) 

(47) 

(48) 


up = 


cos 2 (2iy/|Ao|ct;t) 


(49) 
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For the particular case Aq = 0 we have the following external scale factor 


R(t) oc i 2/(3A ">, 


(50) 


which leads to the following negative def ini te mass density 

KP = ~ iuJV 


(51) 


It is easy to see that, in this particular case A 0 = 0, the external space 
exhibits accelerating expansion provided that 0 < Aw < However, we 
note that it suffers from the negativity of mass density. 

(B) Suppose the gravitational action consists of just the second order 
Euler-Poincare term on its own. This theory is sometimes called the Gauss- 

Bonnet gravity. Let us consider the cases where oo ^ 0. Then the master 

equation becomes 

x = ux 2 : (52) 

which gives the following a power-law expansion 

R(t) oc t 4/(3Aw) , (53) 


and the following mass density 

1 1 

P ~ 2A uAt*' 


(54) 


It is easy to see that in this solution the external space exhibits accelerating 
expansion if 0 < Aw < |. Additionally the mass density takes positive 
values properly provided that A is positive, i.e., A 2 is positive in line with 
string theories. 

Among the explicit solutions we presented above, the case for which the ex¬ 
ternal space expands as R oc | sinh(v4\ak) | 4/,<4iA, E given in (j37|) is the most 
promising solution from the cosmological point of view. We note that in this 
solution the external space behaves as R ~ t 4 /(^ Aw ) for t ~ 0 and R ~ e" 3 A^f * 
for t >> According to this, provided that Aw >3/4, the external space 
evolves from a decelerating expansion phase to the de Sitter expansion as 
the cosmic time t runs. Moreover, if we choose Aw = 2 then the exter¬ 
nal space exhibits exactly the same behavior as the standard ACDM model 

[10] , namely, Aacdm oc sinfE ( \/jtJ, where A is the conventional positive 
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cosmological constant. We note that similar results could be reached in the 
solution (BH) where the second order Euler-Poincare term was omitted. How¬ 
ever, we notice an important difference between these two solutions when we 
consider the corresponding matter densities. The solution (|44|) where the sec¬ 
ond order Euler-Poincare term was omitted suffers from negativity of mass 
density (1461) . while the solution (1HTP where all the Einstein-Hilbert and the 
second order Euler-Poincare terms are considered leads properly to positive 
definite mass density (|38|) . 

5 Concluding Remarks 

We are interested in the idea of a higher dimensional steady state universe 
model because it combines in an interesting way the conventional ideas of a 
Big Bang cosmology with the steady state universe by going to higher di¬ 
mensions. Namely, since the total (3 + n)-dimensional volume of the universe 
is assumed to be constant, the universe in higher dimensions is eternal with 
no beginning and no end. But on the other hand the extra space dimensions 
are expected to contract dynamically to an unobservable size (dynamical 
reduction) while the three dimensional space would be expanding (external 
space), hopefully, in accordance with its evolutionary pattern that we ob¬ 
serve today. The accelerated expansion of the 3-space in this picture would 
be driven by mass that is being continuously transferred from the contract¬ 
ing n-dimensional internal space into the expanding 3-dimensional external 
space. Hence the eternal amount of matter is fixed. It is neither created nor 
exhausted but is redistributed between the external and internal spaces. 

In fact some other higher dimensional cosmological models whose solu¬ 
tions yield a constant higher dimensional volume exist in recent literature 
(See for instance [HI 09], 50] [51] [52]). Yet, the predicted behavior of the 
external space in these papers is not rich enough to explain the observed ac¬ 
celerated expansion, as these solutions exhibit either the power-law expansion 
as in Ref. [16] or de Sitter expansion as in Ref. HZ]. In two recent papers, on 
the other hand, two of the present authors (O.A. and T.D.) considered higher 
dimensional steady state models by imposing the constant volume condition 
from the start and then studied the dynamics implied by this assumption. 
In [46], the authors considered general relativity in (1 + 3 + n)-dimensions 
and obtained various interesting new dynamics for the external space that 
yield a time varying deceleration parameter with oscillating cases included 
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when flat/curved external and curved/flat internal spaces are considered. In 
mi. they extended the study by considering dilaton gravity, that mimics 
the low energy effective string theory action in a very simple manner when 
the dilaton coupling paramater w is fixed to unity. It is found that among 
all possible values in the parameter space of the dilaton coupling parame¬ 
ter and the number of the internal dimensions, the case that corresponds 
to the bosonic string theory (w — 1 and n = 22) gives the best history 
for the effectively four dimensional universe as it fits very well the standard 
history of the universe. Moreover, the case w = 1 with n = 22 not only 
predicts the time and redshift values at which the accelerated expansion has 
started in accordance with the observations, but also predicts the 4 He mass 
fraction in agreement with the mostly used observations for the Big Bang 
Nucleosynthesis. These promising results were obtained using gravitational 
field equations motivated by effective string field theory actions at tree level. 
In this paper within the context of higher dimensional steady state universe 
idea we find it natural to investigate a modified gravity theory motivated 
by loop effects in superstring theories as a quadratic correction at the low 
energy limit. We presented various exact solutions by assuming a higher di¬ 
mensional anisotropic fluid satisfying an equation of state with dimensional 
dichotomy in addition to the (3 + n)-dimensional constant volume assump¬ 
tion. We showed that the external space in our model can exhibit exactly 
the same kinematics with the standard ACDM model even if we omit the 
second order Euler-Poincare term. However the solution where the second 
order Euler-Poincare term was omitted suffers from negativity of mass den¬ 
sity, while the solution where all the Einstein-Hilbert and the second order 
Euler-Poincare terms are considered can maintain the positivity of the mass 
density properly. Hence consideration of the higher order Euler-Poincare 
terms could be crucial for the construction of promising higher dimensional 
steady state universe models. 
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